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ABSTRACT 

The formation of Wolf-Rayet central stars of planetary nebulae ([WR] CSPNe) whose 
spectroscopic appearance mimics massive WR stars remains poorly understood. Least 
understood is the nature and frequency of binary companions to [WR] CSPNe that 
may explain their H-deficiency. We have conducted a systematic radial velocity (RV) 
study of 6 [WR] CSPNe to search for post-common-envelope (post-CE) [WR] binaries. 
We used a cross-correlation method to construct the RV time-series as successfully 
done for massive close binary WR stars. No significant RV variability was detected for 
the late-[WC] type nuclei of Hen 2-113, Hen 3-1333, PMR 2 and Hen 2-99. Significant, 
large-amplitude variability was found in the [WC4] nucleus of NGC 5315. In the [WOl] 
nucleus of NGC 5189 we discovered significant periodic variability that reveals a close 
binary with Porb = 4.04 ±0.1 d. We measured a semi-amplitude of 62.3 ± 1.3 km s“^ 
that gives a companion mass of m 2 ^ 0.5 Mq or m 2 = 0.84 Mq (assuming i = 45°). 
The most plausible companion type is a massive WD as found in Fleming 1. The 
spectacular nebular morphology of NGC 5189 fits the pattern of recently discovered 
post-CE PNe extremely well with its dominant low-ionisation structures (e.g. as in 
NGC 6326) and collimated outflows (e.g. as in Fleming 1). The anomalously long 4.04 
d orbital period is either a once-off (e.g. NGC 2346) or it may indicate there is a 
sizeable population of [WR] binaries with massive WD companions in relatively wide 
orbits, perhaps influenced by interactions with the strong [WR] wind. 

Key words: planetary nebulae: general - stars: Wolf-Rayet - stars: AGB and post- 
AGB - planetary nebulae: individual: NGC 5189 - binaries: spectroscopic 


Planetary nebulae (PNe) are circumstellar gas envelopes 
ejected at the end of the Asymptotic Giant Branch (AGB) 
phase of low-intermediate mass stars of ^IMq to 8Mq. A 
majority of central stars of planetary nebulae (GSPNe) show 
hydrogen rich atmospheres (e.g. Mendez et al. 1988; Mendez 
1991; Napiwotzki & Schoenberner 1995), with their spectra 
containing only weak absorption lines, mainly of hydrogen 


1 INTRODUCTION 


(H) and helium (He). Perhaps the least understood are the 
less numerous emission-line CSPNe with H-deficient atmo¬ 
spheres whose spectral appearance mimics the massive Wolf- 
Rayet (WR) stars. These CSPNe exhibit dense and strong 
stellar winds with mass-loss rates up to ~10“^ M© yr“^ 


* Based on observations made with the South African Astro¬ 
nomical Observatory (SAAO) 1.9-m telescope and the Southern 
African Large Telescope (SALT) under programme 2013-2-RSA- 
005 (PL Miszalski). 


(e.g. Koesterke 2001; Crowther 2008). Their classification is 
similar to massive WR stars where the dominant emission 
lines indicate either a [WC] (He, C and O; Crowther et al. 
1998) or [WN] (He and N; Smith et al. 1996) subtype, where 
the brackets around the spectral type distinguish them from 
their massive counterparts. In contrast to the roughly equal 
division between WN and WC subtypes in massive WR stars 
(van der Hucht 2001), the known examples of [WR] stars are 
heavily skewed towards early and late [WC] types with few 
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intermediate [WC] types (Crowther 2008; see however Gorny 
2014), whereas only two bona-fide [WN] CSPNe were estab¬ 
lished recently in IC 4663 (Miszalski et al. 2012) and A 48 
(Todt et al. 2013; Frew et al. 2014). Comprehending this 
disparity may shed much needed light on how [WR] stars 
form. 

The most peculiar aspect of the newly discovered [WN] 
CSPNe is that their atmospheres are extremely helium rich 
(85-95 per cent by mass, Miszalski et al. 2012; Todt et 
al. 2013). Therefore, they do not fit into the C-rich evo¬ 
lutionary sequence that the vast majority of H-deficient 
post-AGB stars follow (Werner & Herwig 2006). The C-rich 
H-deficient post-AGB evolutionary sequence of [WC] stars 
evolving into PG1159 stars, the latter characterised by sig¬ 
nificantly weaker winds and higher surface gravities, is well 
supported by similar atmospheric abundance patterns be¬ 
tween the two classes (Crowther 2008; Werner et al. 2008). 
Werner (2012) suggested that [WN] CSPNe ‘might’ be the 
progenitors of 0(He) stars, namely a rare group of four high 
surface gravity He-rich post-AGB stars, two of which are 
surrounded by PNe (Rauch et al. 1994, 1996, 1998, 2008; 
Reindl et al. 2013, 2014). However, the existence of [WN] 
CSPNe was still unproven at this time and PB 8 was not 
He-rich enough to be strictly comparable to the 0(He) stars. 
Only with the discovery of the extremely He-rich atmosphere 
of the [WN3] CSPN of IC 4663 (95 per cent He by mass) 
could the 0(He) stars, which are similarly He-rich with He 
mass fractions >90 per cent (Reindl et al. 2013, 2014), be 
securely linked to [WN] CSPNe in a separate He-rich H- 
deficient post-AGB evolutionary sequence [WN]—>-0(He) by 
Miszalski et al. (2012). In this new sequence the 0(He) stars 
are the He-rich analogues of the C-rich PG1159 stars. 

A key question prompted by the discovery of the He-rich 
sequence is whether it operates in a similar fashion to the 
C-rich sequence. Is the appearance of the [WR] phenomenon 
produced by the same recipe? Or is the phenomenon pro¬ 
duced via separate routes with different recipes? There are 
compelling arguments to include binary interactions in the 
formation of some [WC] CSPNe and related objects (De 
Marco & Soker 2002; De Marco 2008) for which the stan¬ 
dard formation scenarios of a late thermal pulse (LTP) (e.g. 
Schoenberner 1979; Blocker 2001; Herwig 2001), a very late 
thermal pulse (VLTP) (e.g. Lawlor & MacDonald 2002), or 
an AGB final thermal pulse (AFTP) (e.g Herwig 2001) may 
not applyQ Miszalski et al. (2012) suggested a merger may 
be responsible for the high He enrichment in IC 4663 and 
Reindl et al. (2013) found a reasonable agreement between 
the atmospheric composition of IC 4663 and that predicted 
from a slow merger of two He-WDs (Zhang & Jeffery 2012). 
A merger scenario would be consistent with the extreme rar- 


^ These scenarios cannot reproduce the moderate level of ~55 per 
cent He by mass in the possible intermediate [WN/WC] CSPN of 
PB 8 (Todt et al. 2010). Miller Bertolami et al. (2011) proposed 
a diffusion induced nova scenario to explain the abundance pat¬ 
tern of PB 8, however these models require low metallicities and 
core masses to reach the extremely He-rich compositions of [WN] 
CSPNe (see the discussion in Todt et al. 2013). While it seems 
that low metallicities may facilitate [WR] formation (Zijlstra et 
al. 2006; Kniazev et al. 2008), the a-element nebular abundances 
of both IC4663 and A 48 are notably Solar, suggesting the Miller 
Bertolami et al. (2011) scenario may not be applicable. 


ity of both [WN] and 0(He) stars compared to their C-rich 
counterparts, however it may not be possible to prove obser- 
vationally. Some other scenarios were discussed by Reindl et 
al. (2014), however we emphasise that many of the scenar¬ 
ios require considerably more theoretical and, in particular, 
observational investigation. 

Unfortunately, there is currently insufficient evidence 
to gauge whether binary interactions in H-deficient post- 
AGB stars are widespread. Are the known binaries merely an 
oddity where the ]WR] component formed independently of 
binary evolution? Or was a binary interaction fundamental 
to forming the ]WR] component? 

Several different possibilities for binary interactions 
were proposed by De Marco & Soker (2002) and De Marco 
(2008), however the paucity of observed binary systems 
makes it difficult to discern to what degree these or some 
other scenarios may play a role. In this work we aim to ad¬ 
dress this by searching for companions in a sample of [WR] 
central stars. It is more practical with modest amounts of 
telescope time to search for close binaries where the orbital 
periods are a few days or less, signifying that a common- 
envelope (GE) interaction has produced the observed short 
orbital periods (Iben & Livio 1993; Ivanova et al. 2013). 
Only two previous post-GE examples have been discov¬ 
ered, namely SDSS J212531.92—010745.9 that consists of a 
PG1159 star with an M dwarf companion (Nagel et al. 2006; 
Schuh et al. 2009) and the [WG7] CSPN of PN G222.8-04.2 
with an undetermined companion type (Hajduk et al. 2010). 
Their literature orbital periods are 0.29 d and 0.63 or 1.26 
d, respectively. 

The majority of close binary CSPNe have been found 
via photometry monitoring including PN G222.8— 04.20 
There are however several limitations with this approach. A 
notable exception is the double degenerate nucleus of Flem¬ 
ing 1 (Boffin et al. 2012), discovered from radial velocity 
(RV) monitoring. The companion in Fleming 1 is a more 
massive WD whose presence is not detectable in the optical 
lightcurve. At least 37 unique early-type [WC], PG 1159 and 
0(He) CSPNe have been subjected to photometric monitor¬ 
ing campaigns to search for non-radial pulsations (Ciardullo 
& Bond 1996; Gonzalez-Perez et al. 2006), however no close 
binaries have been discovered as part of these campaigns. 
If ~15-20 per cent of PNe are expected to host close bi¬ 
naries (Bond 2000; Miszalski et al. 2009a), then we would 
have expected ~5-7 close binaries amongst the sample. It 
is particularly suprising that no binaries were found, espe¬ 
cially since they are amongst the hottest CSPNe known and 
as such they would have produced a very large irradiation 
effect if a close main sequence companion was present (e.g. 
Corradi et al. 2011; Miszalski et al. 2011). If these CSPNe 
were to have companions, they are most likely evolved, mas¬ 
sive WD companions as in Fleming 1 (Boffin et al. 2012). 
Another possibility is that there are main sequence compan¬ 
ions present with orbital periods of several days or months, 
for which an irradiation effect would not be detected (e.g. 
De Marco et al. 2008). An RV monitoring study is therefore 
the least biased method to search for elusive [WR] close bi¬ 
nary CSPNe, however such a study has not been previously 
attempted. Previously researchers have focused primarily on 

2 See Table lAll in Appendix for an updated compilation. 
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Of-type CSPNe, whose photometry (e.g. Handler 2003) and 
line profiles (e.g. Mendez 1990; De Marco et al. 2004; De 
Marco 2009) show significant variability. Perhaps one rea¬ 
son these studies have not been extended to [WR] CSPNe 
is that they were presumed to be similarly afiected by wind 
variability (if not more so for those with faster wind speeds), 
however several studies of massive WR stars have routinely 
found close binaries (Foellmi et al. 2003; Sana et al. 2013; 
Schnurr 2008; Schnurr et al. 2009) and there is no reason 
to suspect the same methods cannot be applied to [WR] 
CSPNe. 

This paper is structured as follows: Section [2] describes 
our sample selection and the observations. Section [3] gives 
an overview of the cross-correlation method used to find RV 
shifts in the spectra. In Section [4] we present our results 
for each object and we perform a ^ analysis to quantify 
variability. Section [5] derives the orbital period and mass 
function of the most variable object, NGC 5189, before dis¬ 
cussing the nature of the companion and the nebular mor¬ 
phology. We conclude in Section |6] 


2 OBSERVATIONS 
2.1 Sample selection 

Table [T] gives the basic properties of the objects considered 
in this pilot study. All were monitored using the Cassegrain 
spectrograph (SpCCD) on the South African Astronomical 
Observatory (SAAO) 1.9-m telescope except for NGC 5189 
which was observed with the Robert Stobie Spectrograph 
(RSS; Burgh et al. 2003; Kobulnicky et al. 2003) on the 
queue-scheduled Southern African Large Telescope (SALT; 
Buckley, Swart & Meiring 2006; O’Donoghue et al. 2006). In 
general we required the spectra to have a resolution of <2A 
(full-width at half maximum, FWHM) and S/N > 40 in the 
continuum to be able to ensure cross-correlation techniques 
could achieve a velocity accuracy of -^10 km s“^ or better. 
Given the relatively poor throughput of SpCCD, we were 
heavily restricted in our choice of [WR] CSPNe to monitor, 
namely those brighter than V < 14 mag. Therefore we se¬ 
lected several suitably bright targets visible at the time of 
observations from Acker & Neiner (2003) and one from Mor¬ 
gan et al. (2001). The magnitude restrictions also imposed 
a bias towards intrinsically brighter late [WC] CSPNe. The 
sample observed by the SAAO 1.9-m is not particularly re¬ 
markable, apart from the peculiar morphologies of Hen 2-113 
(Lagadec et al. 2006) and Hen 3-1333 (Chesneau et al. 2006). 
Both Hen 2-113 and Hen 3-1333 exhibit dual-dust chemistry 
which may have originated from binary evolution. A binary 
system may be necessary to explain the resolved dust disk 
in Hen 3-1333 (De Marco et al. 2002; Chesneau et al. 2006) 
as they are a prevailing feature of post-AGB binaries (Van 
Winckel et al. 2009). The inclusion of NGC 5189 was spe¬ 
cially motivated for SALT observations as described in the 
following. 

NGC 5189 (PN G307.2—03.4) is a bright, remarkably 
peculiar southern PN that has long perplexed astronomers. 
It was discovered in 1826 by J. Dunlop from Australia (Coz¬ 
ens, Walsh & Orchiston 2010), but its unusual appearance 
meant that it was not until much later that it was recognised 
to be a PN with a ‘quasi-planetary’ classification assigned by 


Henize (1967). The filamentary nebula has been described 
as ‘a graceful affair of recurvant gaseous draperies’ (Evans 
1968) and its notorious complexity is the result of multiple 
point-symmetric low-ionisation knots or ansae (Phillips & 
Reay 1983; Hua, Dopita & Martinis 1998; Gonsalves et al. 
2001; Sabin et al. 2012). Phillips & Reay (1983) suggested 
that the point-symmetric features were the result of preces¬ 
sion due to a binary nucleus with an orbital period on the 
order of a few days. Unfortunately, since then no efforts have 
been made in the literature to determine if NGG 5189 has a 
close binary nucleus. The central star has a [WOl] spectral 
type (Crowther et al. 1998) and is amongst the hottest of 
all central stars at 165l_®8 kK (Keller et al. 2014). Ciardullo 
& Bond (1996) monitored the central star for ~9 hours and 
found a low-amplitude pulsation period of 11.51 ± 0.05 min 
that is typical of stars of its ilk (Ciardullo & Bond 1996; 
Gonzalez-Perez et al. 1996), but noted no other variability 
that might otherwise indicate a companion is present. 

2.2 SAAO 1.9-m spectroscopy 

Table [5] gives a summary of the observations with the 
SpCGD spectrograph on the SAAO 1.9-m. They were con¬ 
ducted over two separate weeks namely May 16 to 23 of 2012 
and June 26 to 4 of 2012. We use the modified julian date 
(MJD=JD-2400000.5) as the time of the middle of each ex¬ 
posure. A slit width of 1.5” was used. Both grating 4 (1200 
lines mm“^) and 6 (600 lines mm“^) were used in the ob¬ 
serving runs. A dispersion of 0.49 and 1.09 A/pix for the 
high resolution and low resolution spectra were obtained, 
respectively. The wavelength coverage of grating 4 is 855 A 
at ~ 1.09 A resolution (FWHM) and 1840 A at ~ 2.4 A 
resolution (FWHM) for grating 6. The CCD was read out 
with a binning factor of 1x2 in dispersion and spatial axes, 
respectively. A Copper Argon (CuAr) arc lamp was taken 
before and after each exposure to calibrate the wavelength 
scale. Both short (~ 240s) and long exposures (~ 2400s) 
were taken for each target, since we needed both unsatu¬ 
rated nebular emission lines and deep enough exposures to 
reach a S/N of > 30 - 40 in the continuum. The S/N was 
measured as an average of 4 different S/N at wavelengths ~ 
4220 A, 4755 A, 5290 A and 5620 A for the blue spectra and 
~ 5990 A, 6415 A, 6950 A and 7370 A for the red. The mean 
spectral resolution varied from 2.5 A to 2.9 A . The obser¬ 
vations were mostly made when the seeing ranged between 
~ 1.5” to 2.5”. The wavelength range was chosen according 
to the moon phase (blue wavelengths for grey/dark and red 
wavelengths for bright). 

The 2D frames were trimmed and the overscan region 
bias level was subtracted. The L.A. Cosmic (van Dokkum 
2001) IRAF task was used to cosmic-ray-clean the images. 
Wavelength calibration was performed in the standard fash¬ 
ion and the mean rms achieved in the iraf task identify 
was ~ 0.09 and ~ 0.15 for high and low resolution spec¬ 
tra, respectively. Each solution was applied to the appropri¬ 
ate object spectrum by adding the FITS header keywords 
REFSPECl and REFSPEC2 in the header of the object and 
processing them with the iraf task dispcor. The two arcs 
were interpolated for the final solution, in case there was 
any small shift in the arc during the exposure (e.g. due to 
influence of gravity on the grating causing flexure). We did 
not find evidence of any substantial shift in the arcs. One- 
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Table 1. Basic properties of the sample of [WR] CSPNe observed. 


PNG 

Name 

RA 

Dec. 

V 

Type 

Epochs 

Ref. 

291.3-1-08.4 

PMR 2 

11 34 38.6 

-52 43 33 

13.3 

[WC9/10] 

7 

c 

307.2-03.4 

NGC 5189 

13 33 32.8 

-65 58 27 

14.5 

[WOl] 

14 

b,d 

309.0-04.2 

Hen 2-99 

13 52 30.7 

-66 23 26 

13.3 

[WC9] 

7 

a,b 

309.1-04.3 

NGC 5315 

13 53 57.1 

-66 30 51 

14.4 

[WC4] 

4 

a,b 

321.0-1-03.9 

Hen 2-113 

14 59 53.5 

-54 18 07 

11.9 

[WCIO] 

12 

a,b 

332.9-09.9 

Hen 3-1333 

17 09 00.9 

-56 54 48 

10.9 

[WCIO] 

7 

a,b 


References: (a): Acker and Neiner (2003). (b): Crowther et al. (1998). 

dimensional spectra were extracted using the APALL task. 

The CSPN was extracted, subtracting the immediate nebu¬ 
lar background (4 pixels either side). However, for late [WC] 
types this becomes difficult due to a lack of clearly-defined 
nebular emission lines. 


2.3 SALT RSS spectroscopy 

We obtained 14 RSS spectra of NGC 5189 with the queue- 
scheduled SALT under programme 2013-2-RSA-005 (PI: 
Miszalski). Table [2] includes the log of observations taken 
during January to March 2014. The PG2300 grating was 
used at a camera articulation angle of 70 deg to cover 
A =4442-5462 A which contains several of the highest ion¬ 
isation emission lines from the CSPN of NGC 5189 (e.g. 
O VI, N V, He H) which are formed in the wind closer to 
the [WO] star than lower ionisation species would (e.g. in 
later type [WC] stars), thereby improving the chances of de¬ 
tecting orbital motion (see e.g. Sana et al. 2013). The 1.5 
arcsec wide slit was used at a position angle of 90 deg and 
resulted in a mean resolution of 2.1 A (FWHM) and a mean 
dispersion of 0.32 A/pix. The S/N was measured as an av¬ 
erage of 4 different S/N at wavelengths 4580 A, 4770 A, 
5040 A and 5225 A for all spectra. A ThAr arc lamp expo¬ 
sure was taken before and after each science exposure. Basic 
reductions were applied using the PYSALoj^ package (Craw¬ 
ford et al. 2010) and cosmic ray events were cleaned using 
the LACOSMIC package (van Dokkum 2001). Wavelength cal¬ 
ibration was performed using standard iraf tasks identify, 
REIDENTIFY, FITCOORDS and TRANSFORM by identifying the 
arc lines in each row and applying a geometric transforma¬ 
tion to the data frames. Both the before and after arc expo¬ 
sures were used to derive the fit for the geometric transfor¬ 
mation and the resulting wavelength solution for all spectra 
had an RMS of 0.06 ± 0.02 A. One-dimensional spectra 
were extracted using the APALL task, namely a main extrac¬ 
tion of the CSPN where the immediate nebular background 
(10 pixels on either side corresponding to 5.1 arcsec) was 
subtracted and a second extraction of the nebula were the 
subtracted background was taken from outside the nebula 
(15 pixels on one side only corresponding to 7.6 arcsec) to 
permit measurement of nebular radial velocities close to the 
CSPN position. No flux calibration was performed as the 
spectra were normalised for the cross-correlation analysis 
(see Figure [2]). 

® http://pysalt.salt.ac.za 


(c): Morgan et. al. (2001). (d): Ciardullo et al. (1999). 

In Figures [T] and [2] we show representative spectra of 
each object split over the whole wavelength range. The 
strongest over-subtracted nebular lines in the spectra of 
NGC 5189 have been interpolated over. 

3 METHOD 

3.1 Cross-correlation method &: template 
construction 

The method used for finding radial velocity (RV) shifts in 
the spectra is cross-correlation as described by Foellmi et al. 
(2003). The cross-correlation task used is XCSAO found in 
the RVSAO package (Kurtz & Mink 1998) that is executed in 
the IRAF environment. XCSAO multiplies the Fourier Trans¬ 
form (FT) of the object spectrum with the conjugate of the 
transform of a high S/N ratio zero-velocity template. For a 
detailed overview of the cross-correlation analysis, see Tonry 
& Davis (1979). This method has proved to be very useful in 
finding RV shifts in close massive WR binaries (Foellmi et al. 
2003; Sana et al. 2013; Schnurr 2008; Schnurr et al. 2009). 
Despite the strong winds in [WR] CSPNe, this method is 
still applicable for finding binaries in these systems. 

For a consistent RV monitoring programme, the cross¬ 
correlation template must have high S/N in the continuum. 
The object spectra need to be cross-correlated with this high 
SNR zero-velocity template to obtain accurate RV shifts. 
After having applied a heliocentric correction to the radial 
velocities using the BCVCORR package in iraf, the template 
was constructed. The main steps involved in creating the 
template is outlined as follows: 

(i) Normalise all spectra. 

(ii) Subtract unity from the continuum. 

(iii) Convert all spectra to a logarithmic wavelength scale. 

(iv) Use the spectrum which has the highest S/N in the 
continuum as first template (Tl) for cross-correlation. 

(v) Use the output by XCSAO of the first cross-correlation 
results to shift all spectra to the same “zero-velocity” as the 
Tl template. 

(vi) Finally, combine all the shifted spectra to create a 
zero-velocity template (T2) using the iraf task SCOMBINE 
providing greater weights to the spectra having higher S/N 
in the continuum. 

Once the template T2 was built, the cross-correlation 
is run again on the spectra, but with T2 as the template. 
The spectral regions which contain most of the stellar lines 
(excluding the nebular emission lines as far as possible) were 
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Figure 1. Representative spectra of objects in our sample. 
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Figure 2. Figure 1 (continued). Note the Hell 4686 nebular residual in NGC 5189 and the chip gaps at A = 4955 A to 4961 A and A 
5142 A 5161 A. 
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Table 2. Table showing the log of observations for the sample. 


Object 

MJD 

(Days) 

Date 

(MM/DD/YYYY) 

Exp. time 

(s) 

Grating 

A 

(A) 

Dispersion 

(A pix-l) 

(S/N) 

FWHM 

(A) 

Hen 2-113 

56063.812846 

05/16/2012 

300 

4 

4260-5115 

0.49 

25 

1.5 


56063.823819 

05/16/2012 

900 

?) 

4257-5115 

0.49 

50 

1.1 


56063.839880 

05/16/2012 

900 

J) 

4257-5115 

0.49 

28 

1.1 


56063.869755 

05/16/2012 

1200 

J) 

4257-5115 

0.49 

48 

1.3 


56064.122231 

05/16/2012 

1500 

?) 

4257-5115 

0.49 

31 

1.3 


56064.801191 

05/17/2012 

1500 

?? 

4257-5115 

0.49 

34 

1.3 


56069.996481 

05/22/2012 

1500 

J) 

4257-5117 

0.49 

30 

1.3 


56106.823001 

06/28/2012 

1500 

6 

4025-5936 

1.10 

40 

2.8 


56108.917004 

06/30/2012 

1500 

!) 

4055-5965 

1.10 

80 

2.6 


56108.928125 

06/30/2012 

1500 

?? 

4036-5950 

1.10 

83 

2.7 


56109.813025 

07/01/2012 

1800 

J5 

5722-7565 

1.06 

79 

2.4 


56111.030991 

07/02/2012 

1500 

J5 

5722-7567 

1.06 

53 

2.5 

Hen 2-99 

56068.820887 

05/21/2012 

1800 

6 

4025-5937 

1.10 

43 

2.7 


56106.744197 

06/28/2012 

1800 

J) 

4036-5947 

1.10 

45 

2.7 


56107.745127 

06/29/2012 

1800 

!? 

4036-5950 

1.10 

25 

2.7 


56108.866659 

06/30/2012 

2400 

?) 

5722-7565 

1.06 

29 

2.5 


56109.747499 

07/01/2012 

2400 

!) 

5721-7565 

1.06 

23 

2.5 


56110.976449 

07/02/2012 

2400 

J) 

4034-5948 

1.10 

38 

2.6 


56111.821737 

07/03/2012 

860 

J5 

4033-5947 

1.10 

25 

2.7 

Hen 3-1333 

56106.927875 

06/28/2012 

1200 

6 

4036-5947 

1.10 

61 

2.6 


56108.109317 

06/29/2012 

1800 

!) 

4032-5940 

1.10 

46 

2.7 


56108.899327 

06/30/2012 

1800 

5? 

5722-7567 

1.06 

34 

2.4 


56109.131929 

07/01/2012 

2400 

!? 

5722-7570 

1.06 

50 

2.4 


56109.873401 

07/01/2012 

1500 

J5 

5722-7572 

1.06 

50 

2.4 


56109.996351 

07/01/2012 

1500 

J5 

4035-5948 

1.10 

64 

2.7 


56111.057141 

07/03/2012 

1500 

!? 

4029-5948 

1.10 

89 

2.6 

NGC 5315 

56068.848701 

05/21/2012 

1500 

6 

4025-5939 

1.10 

30 

2.7 


56106.777051 

06/28/2012 

1500 

?) 

4036-5947 

1.10 

24 

2.7 


56109.782759 

07/01/2012 

2400 

JJ 

4033-5944 

1.10 

24 

2.7 


56110.773799 

07/02/2012 

2400 

J5 

4014-5929 

1.10 

23 

2.7 

PMR 2 

56106.715497 

06/28/2012 

1800 

?) 

4036-5948 

1.10 

32 

2.6 


56107.709237 

06/29/2012 

1800 

” 

4036-5948 

1.10 

27 

2.6 


56108.713127 

06/30/2012 

2400 

5? 

5722-7565 

1.06 

40 

2.5 


56109.713509 

07/01/2012 

2400 

5? 

5722-7565 

1.06 

30 

2.4 


56109.843839 

07/01/2012 

2400 

” 

5722-7565 

1.06 

15 

2.6 


56110.721211 

07/02/2012 

1800 

6 

4036-5948 

1.10 

37 

2.7 


56110.735899 

07/02/2012 

2400 


4015-5927 

1.10 

36 

2.7 

NGC 5189 

56672.581101 

01/15/2014 

900 

PG2300 

4440-5462 

0.32 

63 

2.1 


56674.596367 

01/17/2014 

900 

J) 



69 

2.1 


56677.556460 

01/20/2014 

900 

?) 

” 


51 

2.1 


56685.583323 

01/28/2014 

900 

!? 

” 


87 

2.1 


56698.625846 

02/10/2014 

900 

5? 

” 


65 

2.1 


56700.550927 

02/12/2014 

900 

J5 



47 

2.1 


56701.611668 

02/13/2014 

900 

J5 



59 

2.1 


56704.578300 

02/16/2014 

900 

5? 

” 


44 

2.1 


56705.486876 

02/16/2014 

900 

?? 

” 


61 

2.1 


56705.506708 

02/17/2014 

900 

J) 



44 

2.1 


56706.566535 

02/18/2014 

900 

J) 



63 

2.1 


56707.556766 

02/19/2014 

900 

!? 

” 


65 

2.1 


56731.589104 

03/15/2014 

900 

5? 


55 

55 

2.1 


56733.414920 

03/16/2014 

900 

J5 


55 

89 

2.1 


used for cross-correlation because nebular lines may dilute 
any shifts measurable from the stellar lines. Table [3] gives an 
outline of the wavelength ranges used for cross-correlation. 
Following the Foellmi et al. (2003) method we used the bi¬ 
sector method to determine the absolute zero-point velocity 
of T2 to put all the other measurements into heliocentric 


velocities. Table |4] summarises the strongest stellar emission 
lines chosen for this purpose together with their mean he¬ 
liocentric velocity and the mean RV of the nebula measured 
separately from our observations. The online tool vospec 
(Osuna et al. 2005) was used for bisecting the chosen emis¬ 
sion line in the template spectrum and to obtain the absolute 
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Table 3. The wavelength ranges chosen for the cross-correlations 
of stellar lines. The empty fields mean either no observations 
were made in that particular wavelength range or not enough 
spectra were available for cross-correlation. 


Object 

Blue wavelength range 

(A) 

Red wavelength range 

(A) 

Hen 2-113 

4477-4705 & 5050-5500 

6064-6372 & 6736-7050 

Hen 3-1333 

4500-4720 & 4785-5420 

6000-6288 & 6734-7057 

Hen 2-99 

4500-5740 

6630-7300 & 6000-6490 

NGC 5315 

4477-4705 & 5050-5500 

- 

PMR 2 

4477-5400 

6321-6536 & 6838-7055 

NGC 5189 

5200-5450 

- 


RV of the template using the central wavelength (observed) 
and rest wavelength obtained from the NIST websit^. The 
fraction of the emission line height [h) used to obtain the 
middle point varied for different objects. For different lines, 
an estimate was made to the portion which is more or less 
vertical i.e. ~ 0.3/i to O. 8 / 1 . Seemingly, vospec was intro¬ 
ducing a ~0.5 A shift in the template, but we corrected for 
this in our analysis and it does not affect our results. 


3.2 Construction of RV time-series 

The main outputs by XCSAO are the RV shift in km s“^ of the 
object spectrum relative to the T2 template together with 
the height of the cross-correlation peak and the associated 
error in km s“^. In summary, the velocity obtained by the 
bisector represents the mean stellar velocity and the RVs are 
computed using: 


4 RESULTS 
4.1 RV time series 

Tables [BT] and IB2I list the resultant RV time-series for each 
object. Figures |3] and [4] display the data graphically. We 
emphasise that since the [WR] emission lines are formed in 
the wind, the absolute values of the velocities may not be 
meaningful (especially when compared against the nebular 
RVs). To quantify the variability level observed in each ob¬ 
ject, we conducted a analysis of the RV time-series based 
on a null hypothesis which is discussed in detail by Trumpler 
& Weaver (1953). These authors recommend a standard 
computation, from the reduced Xn-i (where n is the num¬ 
ber of RV measurements and n-1 is the number of degrees of 
freedom) which leads to an estimate of the probability that 
the star’s RV is variable. In summary, let’s say the hypoth¬ 
esis to be tested is H : The semi-amplitude “A ” of a radial 
velocity variation equals zero., i.e. there is no statistical dif¬ 
ference between the observed values {xj) and the expected 
values {x). Then alternatively, the hypothesis is that A > 0. 
The test for H is such that: 

(i) Pi is the probability of rejecting H when H is true 
(i.e, saying A > 0, when A = 0). Now, if we set a = 0.05 
arbitrarily, then in this case, Pi ^ a. 

(ii) Pii is the probability of accepting H when H is false 
(i.e, saying A = 0, when A > 0). In this case Pu must be a 
minimum, no matter which orbit. So, maximizing /3 (where 
P = 1 - Pii, is the probability that the test will reject H 
when H is false) in the vicinity of A = 0. 

The compromise consists of allowing the value of P to 
be at least as large as a. Hence, the criterion for rejection 
is: reject the hypothesis that A = 0 whenever: 


Vi — Vxcsao T VSisector + Vh cv 

Where Vxcsao is the RV output by XCSAO, vusector is the 
stellar mean RV calculated using the bisector method and 
Vhcv is the heliocentric velocity correction. 

The error in the individual RVs {vi) was computed in 
quadrature using the equation: 

AUi = y/AVxcsao + ^'^bisector 

Where Avusector is the error measurement for the bi¬ 
sector method and Avxcsao is the error output by XCSAO as 
described by Kurtz & Mink (1998). A peak is selected by 
fitting a smooth curve with accurate values of 5 (the cen¬ 
tral value), h (the height of the central peak) and w (the 
FWHM of the peak). Kurtz & Mink (1998) assumed a sinu¬ 
soidal noise profile, with the halfwidth of the sinusoid equal 
to the halfwidth of the correlation peak. The mean error 
output by XCSAO is a single velocity measurement and is 
given by: 

3 w 

error = - 7- - 7 

8 ( 1 -kr) 

Where w is the FWHM of the correlation peak and r is 
defined in Tonry & Davis (1979), as the the ratio of the 
height of the true peak to the average peak. 


http://www.nist.gov/pml/data/asd.cfm 


(*1 - xf 




2 

a,n — 


1 


where 


• Xa,n-i is the standard value for n — 1 degrees of 
freedom and probability, P = a. 

• a; is the stellar mean (Expected value). 

• Xj is the individual values for the stellar RVs. 

• CTi the individual errors in the RVs. 


The null hypothesis simply states that there is no statis¬ 
tical difference between the observed value (data) and the 
expected value which in our case is the stellar mean. The 
main point is to either accept or reject the null hypothe¬ 
sis with a given probability based on a critical value {a = 
0.05). If the value of the observed is less than Xo.os,n-i! 
we will accept the null hypothesis with a probability of 95%. 
This means that we are 95% sure that the variations in the 
RV are random. Contrarily, if the observed x^ is greater 
than Xo. 05 ,n-i then we will reject the null hypothesis with a 
probability of 95%, meaning the RV is highly variable. The 
convention we will use here is looking at the probability of 
accepting the null hypothesis. 

The above variability test was run on the 6 objects, as¬ 
suming the stellar mean is the expected mean and the results 
are shown in Table [S] For each of them, the Xobs (observed 
Chi-squared) was computed and the probability that the ob¬ 
ject is variable P{Xobs Xn-i) was obtained. For some stars. 
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Table 4. The weighted mean radial velocities for stellar and nebular emission features of each object. The nebular lines used are Ha, 
H/3, H5, H 7 , NII-6583 and OIII-4959. 


Object 

Stellar 

^rest 

(A) 

(RVstellar) 
(km s“^) 

Nebular 

( RVneb) 
(km s“^) 

Hen 2-113 

CII 

4267.26 

-60±4 

Ho, H/3, H5, H 7 , Nil 

-63±2 

Hen 3-1333 

CHI 

5695.19 

-43±4 

Ho, H/3, H 7 , Nil 

- 66±10 

Hen 2-99 

CHI 

5695.19 

-56±12 

Ho, H/3, Nil 

-95±13 

NGC 5315 

Hel 

4471.47 

84±15 

Ho, H/3, H 7 , Nil, OIII 

-24±8 

PMR 2 

CII 

4267.18 

4±5 

Ho, H/3, Nil 

-66±23 

NGC 5189 

OVI 

5290.65 

23±40 

Hp 

-8±4 






Hen 3-1333 


1 — p — 

1 NMl 1 

.— ^ 

- 



107.0 107.5 108.0 108.5 109.0 



MJD-56000 (Days) 


Figure 3. Measured RV time-series for the nuclei of Hen 2-113, Hen 2-99 and Hen 3-1333. 
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Figure 4. Measured RV time-series for the nuclei of NGC 5315, PMR 2 and NGC 5189. 


the RV variability was so large that it becomes obvious with 
only a few observations to decide whether it is a binary or 
not. However, in most cases we lack sufficient observations 
to properly constrain the amount of variability in the RV. 
In cases where small random fluctuations in the RV from 
the mean occur, which we would expect due to instrumental 
uncertainties or wind variability, we would identify the RV 
as being constant. However, a highly variable RV from the 
mean would be one having large enough RV variations that 
cannot be assigned solely to random variations. In the case 


of a close binary CSPN we require a statistically significant 
periodic signature in the RV time-series to be present. 

The objects observed with the SAAO 1.9-m are not 
suffiently variable to prove a binary companion is present 
at this time, however two objects stand out from the rest. 
Hen 2-99 has stellar RV shifts with a peak-to-peak value of 
the order 30 km s“^ from the mean and a hint of periodicity 
is seen with a peak-to-peak of ~ 53 km s~^. A sinusoid fits 
reasonably well with the RV of Hen 2-99 with a period of 
5.3 d. More data is needed to confirm this period. Similarly, 
NGC 5315 shows much higher RV variability (Pnuii = 1 %), 
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Table 5. Results of the variability test carried on individual ob¬ 
jects. Column 1 shows the object name, the 2”'^ column is the 
number of RV measurements (degrees of freedom+1), the S’’'* 
column is the mean RV, the 4*^ column is the standard deviation 
in the RV, the 5*^ column is the observed s-nd column 6 shows 
the probability of accepting the null hypothesis based on the 
value. 


Object 

n 

(RV) 
(km s ^) 

(T 

(km s“^) 

Xobs 

Pnull 

Hen2-113 

18 

-59.5 

4.1 

8.8 

94% 

Hen3-1333 

15 

-44.1 

4.4 

0.2 

100 % 

Hen2-99 

8 

-55.7 

12.0 

4.5 

72% 

NGC 5315 

4 

87.3 

14.0 

15.3 

1 % 

PMR 2 

7 

4.0 

4.8 

3.0 

80% 

NGC 5189 

14 

23.2 

40.1 

690 

0 % 


but we have only 4 epochs for this object. The non-periodic 
and relatively small amplitude RV shifts in the remaining 
objects that we see are most probably due to stochastic wind 
variability (e.g. Grosdidier et al. 2000, 2001). 

NGC 5189 is the only object in the sample that shows 
very high RV variability with a peak-to-peak velocity of the 
order of 120 km s (see Fig. 3)). From the variability 
analysis in Table [5] and Fig. |4] it is quite clear that it is 
the only object which shows 100% variability. The following 
section further analyses this unique dataset. 


5 THE CLOSE BINARY CS OF NGC 5189 


The second panel in Figure [5] is a sinusoid model deter¬ 
mined from fitting the data folded on the 4.04 day period in 
phase space. We used the mpfit and mpfitfun idl codes de¬ 
signed to perform a Levenberg-Marquardt least-squares fit 
of a user supplied model to a function (Markwardt 2012). 
The main fitting was done by mpfit after specifying a set 
of initial parameters. Given the data and their uncertain¬ 
ties, MPFITFUN computes the best set of model parameters 
which match the data and returns them in an array. In time 
space we determined a 4.04 day period fit of the form A 
sin(27r/ — (fi) + B to the data where A denotes the ampli¬ 
tude, f is the frequency fixed at 0.247 day“^, cf) is the phase 
and B is the mean, with A = 62.3 km s“^, (j) = 0.43 and 
B — 14.0 km s“^. This corresponds to an ephemeris deter¬ 
mined from the minimum radial velocity at MJD (min RV) 
= (2456700.24 ± 0.02) -b (4.04 ± 0A)E. The error in the 
orbital period was determined by fitting a gaussian with a 
mean of 4.04 d and a standard deviation of 0.1 d to the peak 
period in the periodogram. The bottom two panels of Fig. 
[5] shows this fit along with the RV time-series data showing 
the reasonable agreement to our observations. Figure |S] gives 
another graphical representation of the periodic motion in 
the RV time-series data of NGC 5189. The trail diagrams are 
plotted over 2 phases using the motion of the OVI-5290 stel¬ 
lar emission line and, for comparison, the static nebular 
emission line. 

In the case of close-binaries with a good sinusoidal fit, it 
is fair to assume a low (e<Cl) or zero eccentricity. In this case 
the mass function of two stars of masses mi and m 2 , where 
only one of the stars is visible in the spectrum (single-line 
spectroscopic binaries) is described by: 


5.1 Orbital period analysis and mass function 

Since our data is not evenly sampled, we used the Lomb- 
Scargle (Lomb 1976; Scargle 1982) method to search for an 
orbital period. The Lomb-Scargle method is mainly based 
on a least-squares fit of sinusoids (Press & Rybicki 1989). 
Figure [ 5 ] shows four significant peaks in the periodogram at 
periods 4.04 days, 3.54 days, 1.32 days and 0.80 days. How¬ 
ever, we find that the 3.54 day, 1.32 day and 0.80 day peaks 
are most likely to be aliases, as they do not persist when we 
ran the Lomb-Scargle for a second time after extracting the 
4.04 day period. Moreover, to test the validity of the peri¬ 
ods, we further analysed the data according to the method 
described by Tanner (1948). Assuming the actual period is 
4.04 d, both the 1.32 d and the 0.8 d periods satisfy Equa¬ 
tion 2 in Tanner (1948). Based on the relatively high value 
of 8.3 for the 3.54 day period sinusoid fit, we further re¬ 
ject this period. This leaves us with the best period being 
4.04 d which has the lowest x^ of 3.0 out of all periods. 

A significance test was also carried out on the peri¬ 
ods found using a Monte-Carlo simulation, where a random 
number generator was used to scramble the original radial 
velocities within one standard deviation, creating 10000 ran¬ 
domised RV curves of the same structure. A Lomb-Scargle 
analysis was then performed on the randomised RV curves 
and the highest powers were recorded for each of them. The 
powers were then sorted in ascending order after which the 
90%, 99% and 99.9% significance levels were obtained (dis¬ 
played as horizontal lines in Fig. 0. The only period which 
was found to be above the 99.9% level and acceptable as per 
the sinusoid model fit was the 4.04 day. 


/(mi, m 2 ) 


ml 


(mi -I- m2)2 



where P is the orbital period (4.04 ± 0.1 d), i is the incli¬ 
nation of the orbital plane to the line of sight, and v\ is the 
semi-amplitude of the RV curve (ui = 62.3 ± 1.3 km s“^). 

Table [6] lists companion masses m 2 calculated for a 
range of inclination values and an adopted mi mass of 0.596 
Mq for the [WOl] primary. The value of mi was calcu¬ 
lated as an average of the masses of Miller Bertolami & Al- 
thaus (2006) evolutionary tracks either side of the location 
of NGC 5189 in figure 13 of Keller et al. (2014). There are 
several caveats associated with this assumption. Firstly, the 
mass estimates derived from interpolating between evolu¬ 
tionary tracks are distance dependent and distances to PNe 
are notoriously uncertain. Keller et al. (2014) adopted a dis¬ 
tance of 0.55 kpc to NGC 5189 (Stanghellini et al. 2008), 
but the distance may be as large as 1.44±0.27 kpc based on 
extinction and kinematic distance estimates (Frew 2008). 
The luminosity would therefore increase from ~ 2.73 Lq 
(K eller et al. 2014) to ~ 3.86 Lq (Frew 2008), resulting in 
a mass closer to -^0.62 Mq. However, it is not necessarily a 
matter of adopting this mass as this method of estimating 
CSPNe masses does not agree well with spectroscopic mass 
estimates for CSPNe (see section 4.2.3 of Moe & De Marco 
2006 and ref. therein), underlying the intrinsic unreliability 
of such model-dependent methods. Finally, the evolutionary 
tracks assume single star evolution which may not apply to 
a post-CE binary such as the nucleus of NGC 5189. 

In summary, these caveats make it difficult to choose 
any ‘best’ value for mi. As the orbital inclination of the neb- 
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MJD-56000 (Days) 


Figure 5. Top: The Lomb-Scargle periodogram showing the significance levels of 90%, 99% and 99.9%. Middle: The radial velocity 
folded on the 4.04 day period fitted with a sinusoidal model. Bottom: Plot of the radial velocities folded on the 4.04 day period fitted 
with a sinusoid. 
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Figure 6. Left: Trail diagram of the nebular H/3 line. Right: Trail diagram of the stellar OVI-5290 emission line. The dashed line shows 
the sinusoidal fit to the RV curve in phase space. 


ula is currently unconstrained (Sabin et al. 2012), it is only 
possible in this work to give an indication of plausible com¬ 
panion masses. The unknown inclination introduces a much 
larger uncertainty in m 2 than the small uncertainty in mi. 
Therefore, for the purposes of this work, it is reasonable to 
retain the adopted value of mi = 0.596 Mq that is essen¬ 
tially the same as the standard mass assumed for CSPNe 
of 0.6 Mq (e.g. Crowther et al. 2006; Crowther 2008). Even 
if a slightly more massive mi were adopted, it would only 
result in a small increase in m 2 . 


5.2 Nature of the companion 

We give a summary of the properties of the CS of NGC 5189 
in Tableland Table [8] gives the orbital parameters derived 
in this work. The mass function puts a lower limit on m 2 (0.5 
Mq ± 0 .01) which is reached only if the system is edge-on 
with i = 90°. For inclinations less than 40 degrees we hnd 
improbably high values of m 2 , suggesting these inclinations 
can be discarded. At other inclinations the possible compan¬ 
ions are a main sequence star or a WD. It may be possible 
for an orbital period of 4.04 d that a main sequence compan¬ 
ion could produce an irradiation effect, especially given the 
165 kK temperature of the [WOl] component. This effect 


Table 6. Companion masses (m 2 ) for NGC 5189 for a variety of 
orbital inclination angles. 


Inclination 

(°) 

7712 

(Mq) 

Am2 

(Mo) 

30 

1.56 

0.03 

35 

1.21 

0.02 

40 

0.99 

0.02 

45 

0.84 

0.01 

50 

0.74 

0.01 

55 

0.67 

0.01 

60 

0.61 

0.01 

65 

0.57 

0.01 

70 

0.54 

0.01 

75 

0.52 

0.01 

80 

0.51 

0.01 

85 

0.50 

0.01 

90 

0.50 

0.01 


could be up to 0.1-0.2 mag (De Marco et al. 2008), but this 
has not yet been observed (e.g. Ciardullo & Bond 1996). 
There is also no evidence to suggest a cool companion is 
present since the intrinsic {V — I) colour of —0.32 ± 0.08 
is typical of hot blue CSPNe {(V — I)o ^ —0.4, Ciardullo 
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Table 7. Properties of NGC 5189 and the [WOl] primary. 


PN G 

307.2-03.4 


c{Hfl) 

0.47 ± 0.08 

Garcia-Rojas et al. (2012) 

V (mag) 

14.53 

Ciardullo et al. (1999) 

I (mag) 

14.35: 

Ciardullo et al. (1999) 

1 

o 

-0.32+0.08 

This work 

Type 

[WOl] 

Crowther et al. (1998) 

Voo (km s“l) 

2500+250 

Keller et al. (2014) 

Teff (kK) 

165t« 

Keller et al. (2014) 

mi (Mq) 

0.596 

Keller et al. (2014) 


Table 8 . Orbital parameters of NGC 5189. 


Thrb(^) 

4.04 + 0.1 

e 

0 (fixed) 

K (km s“^) 

62.3 + 1.3 

7 (km s“^) 

14 + 1 

m2 (Mq) 

^ 0.5 

To (MJD) 

56699.74+0.02 


et al. 1999). The intrinsic faintness of the primary, namely 
Mv = +2.7 to +4.8 mag for d = 1.44 kpc (Frew 2008) or 
0.55 kpc (Stanghellini et al. 2008), respectively, also means 
that it is unlikely that a cool companion is hiding in the 
glare of the primary. Although the orbital inclination is cur¬ 
rently unconstrained (Sabin et al. 2012), we can estimate 
an intermediate inclination of 45 degrees based on the ap¬ 
parent nebular morphology, i.e. the relatively symmetrical 
inner nebula with a pinched waist (see Sect. 15.31) 1 and that 
the collimated outflows are neither completely inside or out¬ 
side the nebula. The companion is therefore likely to be a 
slightly more massive WD of ~0.9 Mq, comparable to the 
companion discovered in Fleming 1 (Boffin et al. 2012). 

The 4.04 d orbital period of NGC 5189 is one of the 
longest currently known amongst post-CE CSPNe (Fig. [7]) 
and falls into the zone with periods >1 day where there is a 
substantial deficit of post-CE binaries compared to CE pop¬ 
ulation synthesis models (Rebassa-Mansergas et al. 2008; 
Miszalski et al. 2009a; Davis et al. 2010; Nebot Gomez- 
Moran et al. 2011). There are two possible interpretations 
for the position of NGC 5189 in the post-CE orbital period 
distribution. One is that the longer period is a one-off co¬ 
incidence as in NGC 2346 and this would be in line with 
the observed deficit. The other possibility is that NGC 5189 
reflects the tip of an iceberg where most [WR] close bina¬ 
ries would have similarly long orbital periods (or longer). 
This could be explained by the stronger and more extended 
wind of the [WR] component that could interact with the 
companion (e.g. via wind Roche-lobe overflow, Mohamed & 
Podsiadlowski 2007, 2011). If the interaction were to facili¬ 
tate the CE phase at greater orbital separations than in the 
average post-CE CSPN (WD and main sequence), then the 
resultant orbital period might be longer than the average. 
At present with only two [WR] post-CE CSPNe known, it is 
difficult to decide which of these two possibilities is correct. 



Figure 7. Orbital period distribution of published close binary 
GSPNe (see Appendix ETJ. 

5.3 The post-CE nebular morphology 

The main reason for including NGC 5189 in our study was 
that it contains all the hallmarks of post-CE nebulae as first 
outlined by Miszalski et al. (2009b). These hallmarks in¬ 
clude low-ionisation structures (LIS; see e.g. Gongalves et 
al. 2001), collimated outflows or jets and bipolar nebulae 
(see also Miszalski et al. 2011b and Miszalski 2012). Ad¬ 
ditional recent discoveries of post-CE nuclei in other PNe 
have further reinforced these hallmarks (Corradi et al. 2011; 
Miszalski et al. 2011a,b,c; Boffin et al. 2012; Jones et al. 
2014). The abundance of LIS filaments in NGC 5189 (e.g. 
Sabin et al. 2012), the majority of which point to the cen¬ 
tral star, is exactly what we see in NGC 6326 (Miszalski 
et al. 2011b) and similarly high levels of LIS are also seen 
in NGC 6778 (Miszalski et al. 2011b; Guerrero & Miranda 
2012) and Hen 2-11 (Jones et al. 2014). Secondly, Sabin et al. 
(2012) demonstrated that the outermost LIS of NGC 5189 
have the most extreme velocities of the nebula. Combined 
with the central S'-shaped feature, these outermost LIS are 
symptomatic of a processing outflow (e.g. Lopez et al. 1993) 
that is best explained by jets being launched from a precess- 
ing accretion disk around a companion (Cliffe et al. 1995; 
Raga et al. 2009; Boffin et al. 2012). Furthermore, leaving 
the LIS aside, the intrinsic morphology of NGC 5189 appears 
to be a bipolar outflow best seen in the [O III] emission line. 
It is pinched at the waist at the edges of the minor axis that 
is closely aligned with slit position 4 in figure 2 of Sabin et 
al. (2012). As emphasised in Miszalski et al. (2011b) with 
NGC 6326, there is a stark disconnect between the intrinsic 
morphology of the nebula and the LIS, the latter of which 
may form separately to the diffuse components. 


6 CONCLUSIONS 

We carried out an RV monitoring program of 6 Galac¬ 
tic [WR] CSPNe using a well established cross-correlation 
method successfully used in massive WR star studies 
(Foellmi et al. 2003). Four of the CSPNe in our sample are 
late-types ([WCL[) namely Hen 2-99, Hen 2-113, Hen 3-1333 
and PMR 2, and two are early-type ([WCE]), namely the 
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CSPNe of NGC 5189 and NGC 5315. Our main conclusions 
are as follows: 

(i) No significant variability was detected in Hen 2-113, 
Hen 3-1333 and PMR 2. Hen 2-99 may be variable with 
a putative periodicity of 5.3 d, but more observations are 
needed to further investigate. Similarly, NGC 5315 shows a 
much greater level of variability but will again require more 
observations to search for an orbital period. 

(ii) NGC 5189 showed significant RV variability in stark 
contrast to the other objects in our sample. A significant 
period of 4.04 d was found from the Lomb-Scargle peri- 
odogram and fitting the phased data gave a peak-to-peak 
amplitude of ~ 124 km s“^. The orbital motion is clearly 
seen in the trailed diagram of the stellar O VI 5290 emis¬ 
sion line compared against the static H/3 nebular emission 
line. The mass of the companion m2 has a lower limit of 0.5 
M0 and may be ~0.9 Mq for a reasonable assumption of 
the nebular inclination of i = 45 degrees. A more massive 
WD companion is likely to be present, explaining the lack of 
an irradiation effect in previous photometric observations of 
the object. Further spatio-kinematic modelling of NGC 5189 
as in Sabin et al. (2012) is strongly encouraged to further 
constrain the inclination and therefore companion mass. 

(iii) The spectacular nebula of NGC 5189 further fits the 
trend for post-CE nebulae to be dominated by low-ionisation 
structures (e.g. NGC 6326, Miszalski et al. 2011) and to pos¬ 
sess precessing outflows (e.g. Fleming 1, Boffin et al. 2012) as 
first outlined by Miszalski et al. (2009b). Several of the PNe 
with [WR] nuclei discussed in Miszalski et al. (2009b) are 
therefore excellent candidates for RV monitoring to discover 
new [WR] binaries. 

(iv) The discovery of a second close binary system con¬ 
taining a [WR] component strongly suggests mergers are not 
involved in the formation of most [WC[ CSPNe. It remains 
to be seen whether [WN[ CSPNe are binary systems. Indeed, 
the relatively long orbital period of NGC 5189 could be ei¬ 
ther a one-off coincidence (e.g. NGC 2346) or alternatively 
it could indicate that potentially many more [WR] bina¬ 
ries may be found if appropriate RV monitoring surveys are 
conducted. We speculate that wind interactions between the 
[WR] component and its companion (e.g. wind Roche-lobe 
overflow) may be responsible for the longer period. Further 
RV monitoring of [WR] CSPNe is urged to try and discover 
more systems to place constraints on the formation of [WR[ 
CSPNe. 
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APPENDIX A: LIST OF CLOSE-BINARIES 

Table ED gives a list of published post-CE PNe updating 
Miszalski et al. (2011c). In the Miszalski et al. (2009a) sam¬ 
ple, there are still some objects that require spectroscopic 
confirmation that the variable star identified is the CSPN. 
Miszalski (2009) obtained Gemini CMOS spectroscopy for 
several in the sample. Some objects were removed based 
partially on the spectra (Miszalski et al. 2011c) and several 
other objects were confirmed by Miszalski (2009) as true cen¬ 
tral stars (K 6-34, H 2-29, BMP 1800-3408, PPA 1759-2834, 
Pe 1-9, PPA 1747-3435, PER 1757-2824, PER 1756-3342, 
Sab 41, M 2-19 and M 3-16). Unpublished spectra of other 
objects have also been obtained, but these will be discussed 
elsewhere. We have removed Te 11 which is a cataclysmic 
variable (Drake et al. 2014) and MPA 1508-6455 which re¬ 
quires further work to prove its long 12.5 d period. A period 
of 1.26 d is adopted for PN G222.8—04.2 based on unpub¬ 
lished SALT RSS radial velocities from our 2013-2-RSA-005 
programme. 
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Table Al. An updated list of 42 close binary CSPNe whose status is fairly well established. 


PN G 

Name 

Period Discovery reference 
(days) 

053.8-03.0 

Abell 63 

0.46 

Bond et al. 1978 

215.6+03.6 

NGC 2346 

15.99 

Mendez & Niemela 1981 

009.6+10.5 

Abell 41 

0.23 

Grauer &C Bond 1983 

055.4+16.0 

Abell 46 

0.47 

Bond 1985 

283.9+09.7 

DS 1 

0.36 

Drilling 1985 

136.3+05.5 

HFG 1 

0.58 

Grauer et al. 1987 

253.5+10.7 

K 1-2 

0.68 

Bond Grauer 1987 

005.1-08.9 

Hf 2-2 

0.40 

Lutz et al. 2010 

017.3-21.9 

Abell 65 

1.00 

Bond Livio 1990 

329.0+01.9 

Sp 1 

2.91 

Bond Livio 1990 

355.2-03.6 

HaTr 4 

1.74 

Bond Livio 1990 

144.8+65.8 

BE UMa 

2.29 

Liebert et al. 1995 

135.9+55.9 

SBS 1150+599A 

0.16 

Tovmassian et al. 2004 

341.6+13.7 

NGC 6026 

0.53 

Hillwig et al. 2010 

349.3-01.1 

NGC 6337 

0.17 

Hillwig et al. 2010 

359.1-02.3 

M 3-16 

0.57 

Miszalski et al. 2008 

357.6-03.3 

H 2-29 

0.24 

Miszalski et al. 2008 

000.2-01.9 

M 2-19 

0.67 

Miszalski et al. 2008 

005.0+03.0 

Pe 1-9 

0.14 

Miszalski et al. 2009a 

355.3-03.2 

PPA 1747-3435 

0.22 

Miszalski et al. 2009a 

355.7-03.0 

H 1-33 

1.13 

Miszalski et al. 2009a 

354.5-03.9 

Sab 41 

0.30 

Miszalski et al. 2009a 

000.6-01.3 

B1 3-15 

0.27 

Miszalski et al. 2009a 

359.5-01.2 

JaSt 66 

0.27 

Miszalski et al. 2009a 

358.7-03.0 

K 6-34 

0.20 

Miszalski et al. 2009a 

357.0-04.4 

PHR 1756-3342 

0.26 

Miszalski et al. 2009a 

001.8-02.0 

PHR 1757-2724 

0.80 

Miszalski et al. 2009a 

001.2-02.6 

PHR 1759-2915 

1.10 

Miszalski et al. 2009a 

005.0-03.1a 

MPA 1759-3007 

0.50 

Miszalski et al. 2009a 

001.9-02.5 

PPA 1759-2834 

0.31 

Miszalski et al. 2009a 

357.1-05.3 

BMP 1800-3408 

0.14 

Miszalski et al. 2009a 

000.9-03.3 

PHR 1801-2947 

0.32 

Miszalski et al. 2009a 

222.8-04.2 

PHR 0654-1045 

1.26 

Hajduk et al. 2010 

054.2-03.4 

The Necklace 

1.16 

Gorradi et al. 2011 

068.1+11.0 

ETHOS 1 

0.53 

Miszalski et al. 2011a 

049.4+02.4 

Hen 2-428 

0.18 

Santander-Garcia et al. 2015 

034.5-06.7 

NGC 6778 

0.15 

Miszalski et al. 2011b 

338.1-08.3 

NGC 6326 

0.37 

Miszalski et al. 2011b 

290.5+07.9 

Fleming 1 

1.19 

Boffin et al. 2012 

259.1+00.9 

Hen 2-11 

0.61 

Jones et al. 2014 

086.9-03.4 

On 5 

0.36 

Gorradi et al. 2014 

307.2-03.4 

NGC 5189 

4.04 

This work 
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APPENDIX B: LOG OF 
CROSS-CORRELATION RESULTS 
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Table Bl. Log of cross-correlation wavelength ranges, RV shifts, errors and XCSAO cross-correlation peak heights h. 


Object 

A range 

(A) 

Grating 

MJD (mid) 
(Days) 

Vi 

(km s“^) 

Avi 

(km s“^) 

h 

Hen2-113 

4477-4705 

4 

56063.812846 

-59.7 

4.9 

0.97 




56063.823819 

-60.0 

4.8 

0.97 




56063.836879 

-60.1 

4.7 

0.97 




56063.869755 

-60.3 

4.5 

0.98 




56064.122231 

-56.0 

4.0 

0.97 




56064.801191 

-59.0 

5.0 

0.97 


5050-5500 


56063.812846 

-61.4 

7.0 

0.81 




56063.823819 

-60.0 

6.0 

0.82 




56063.836879 

-63.0 

6.0 

0.89 




56063.869755 

-60.0 

6.0 

0.92 




56064.122231 

-61.1 

7.2 

0.90 




56064.801191 

-52.9 

7.0 

0.74 


5050-5500 

6 

56069.996481 

-50.5 

5.0 

0.97 




56106.823001 

-60.1 

4.9 

0.98 




56108.917004 

-64.9 

6.6 

0.92 




56108.928125 

-58.2 

5.6 

0.91 




56109.813025 

-68.5 

5.7 

0.89 




56111.030991 

-56.2 

4.9 

0.98 

Hen3-1333 

4500-4720 

6 

56106.927875 

-41.8 

4.4 

0.98 




56108.109317 

-40.2 

4.4 

0.98 




56109.996351 

-48.2 

4.7 

0.98 




56111.057141 

-47.3 

4.3 

0.99 


4785-5420 


56106.927875 

-44.8 

4.1 

0.98 




56108.098274 

-47.4 

13.0 

0.58 




56108.109317 

-35.0 

4.3 

0.98 




56109.996351 

-47.1 

4.1 

0.98 




56111.057141 

-49.1 

4.6 

0.97 


6000-6288 


56108.899327 

-45.0 

4.0 

0.98 




56109.131929 

-35.0 

5.3 

0.98 




56109.873401 

-45.1 

3.9 

0.98 


6734-7057 


56108.899327 

-40.5 

3.9 

0.98 




56109.131929 

-41.2 

4.3 

0.98 




56109.873401 

-45.1 

3.9 

0.98 

Hen2-99 

4500-5740 

6 

56107.745127 

-35.6 

13.7 

0.98 




56068.820887 

-59.0 

12.7 

0.98 




56110.976449 

-56.0 

14.3 

0.98 




56111.821737 

-54.5 

12.6 

0.98 




56106.744197 

-44.5 

15.0 

0.97 


6630-7300 


56108.866659 

-62.0 

14.1 

0.96 




56109.747499 

-75.5 

16.2 

0.96 


6000-6490 


56108.866659 

-58.6 

15.0 

0.87 
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Table B2. Table [ bD (continued). 


Object 

A range 

(A) 

Grating 

MJD (mid) 
(Days) 

Vi 

(km s“^) 

Avi 

(km s“^) 

h 

NGC 5315 

4477-4705 

6 

56068.848701 

93.9 

6.5 

0.76 




56106.777051 

95.8 

11.5 

0.83 




56110.773799 

65.8 

6.1 

0.82 


5050-5500 


56109.782759 

79.7 

6.4 

0.96 

PMR 2 

4477-5400 

6 

56107.709237 

2.1 

6.4 

0.98 




56110.721211 

4.6 

12.0 

0.80 




56110.735899 

-1.3 

6.2 

0.98 




56106.715497 

0.4 

6.1 

0.98 


6321-6536 


56108.713129 

9.5 

6.4 

0.96 




56109.713509 

9.7 

8.5 

0.94 




56109.843839 

10.0 

8.4 

0.84 

NGC 5189 

5169-5443 

PG 2300 

56672.581101 

-34.6 

6.0 

0.90 




56674.596367 

36.1 

7.9 

0.85 




56677.556460 

56.7 

6.0 

0.86 




56685.583323 

61.6 

5.9 

0.85 




56698.625846 

70.3 

5.5 

0.85 




56700.550927 

-40.4 

5.6 

0.87 




56701.611668 

51.6 

4.4 

0.90 




56704.578300 

-37.0 

6.8 

0.88 




56705.486876 

49.4 

3.7 

0.88 




56705.506708 

44.1 

6.1 

0.83 




56706.566535 

57.6 

5.1 

0.87 




56707.556766 

-14.9 

4.9 

0.91 




56731.589104 

20.8 

4.7 

0.85 




56733.414920 

4.2 

3.7 

0.90 







